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We study the effect of uniform uniaxial strain on the ground state electronic configuration of 
a thin film manganite. Our model Hamiltonian includes the double-exchange, the Jahn- Teller 
electron-lattice coupling, and the antiferromagnetic superexchange. The strain arises due to the 
lattice mismatch between an insulating substrate and a manganite which produces a tetragonal 
distortion. This is included in the model via a modification of the hopping amplitude and the 
introduction of an energy splitting between the Mn e 9 levels. We analyze the bulk properties of 
half-doped manganites and the electronic reconstruction at the interface between a ferromagnetic 
and metallic manganite and the insulating substrate. The strain drives an orbital selection modifying 
the electronic properties and the magnetic ordering of manganites and their interfaces. 

PACS numbers: 75.47.Gk, 75.10.-b, 75.30.Kz, 75.50.Ee 



I. INTRODUCTION 

Manganites are strongly correlated oxides that show 
a large variety of magnetic and electronic phases due 
to a strong interrelation be twee n the orbital, charge 
and spin degrees of freedom!^ They are particularly 
well known for the measured colossal magnetoresistance 
and the half-metallicity, which is responsible for the 
large tunneling magnetoresistance observed in mangan- 
ite/insulator/manganite trilayers!^ The current inter- 
est on oxide heterostructures^ and the electronic recon- 
struction occurring at their interfaces"^ has also been re- 
flected on a renewed interest on heterostructures involv- 
ing manganites with different properties. 10 " 22 In differ- 
ent all-manganite heterostructures, it has been observed 
that the properties of thin manganite layers may be mod- 
ified with respect to their bulk behavior, for instance, by 
the appearance of a ferromagneti c mome nt in a nomi- 
nally antiferromagnetic manganit e ! 15 * 16 * 22 ^ or by the for- 
mation of a ferromagnetic two dimensional electron gas 
at the interfac e betwee n two antiferromagnetic insulating 
manganites ESHUl 14 l 1 7tf H] Orbital reconstruction, a mod- 
ification of the orbital occupancy at interfaces between 
different materials, has also been observedP2HH] 

Manganites have the pseudocubic perovskite structure 
with chemical composition Ai-^A'^MnOa with A typ- 
ically a trivalent rare-earth (e.g. La or Pr) and A' a 
divalent cation (e.g. Ca or Sr). 1 — x is the concentration 
of electrons moving on the Mn e g (x 2 — y 2 and 3z 2 — r 2 ) 
orbital bands. Mn ions are in the center of oxygen oc- 
tahedra that may undergo Jahn- Teller (JT) distortions. 
These lattice distortions couple to the charge and orbital 
degrees of freedom producing a splitting of the e g levels 
sometimes associated with charge and/or orbital order- 
ing. The competition between the JT, the antiferromag- 
netic (AF) superexchange, and the kinetic energy via the 
double exchange (DE) interaction (Hund's coupling is as- 
sumed to be infinite), leads to a complex phase diagram 
as a function of composition and doping! 2 ^ In particular, 
at half-doping (x = 0.5) many manganites are insulating 
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FIG. 1: (Color online) Cartoon of the possible spin and orbital 
order configurations for half-doped manganites. FM labels 
a three dimensional ferromagnet while the other three (A, 
C, and CE) are different antiferromagnetic orderings. The 
FM configuration is orbital disordered (represented here by 
isotropic spherical orbitals). The A-type, with FM planes 
coupled antiferromagnetically, favors the occupancy of the 
x 2 — y 2 orbitals. The C-type, with FM lines in the z-direction 
coupled antiferromagnetically, favors the occupancy of the 
3z 2 — r 2 orbitals. The CE-type ordering consists of FM zig- 
zag chains coupled antiferromagnetically (only the xy plane is 
shown here) and is associated with a peculiar orbital ordering 
(alternating 3a; 2 — r 2 , 3j/ 2 — r 2 and more isotropic orbitals) 
and checkerboard charge ordering. 



and show charge, orbital and antiferromagnetic ordering 
(of the CE-type, see Fig. [TJ) while for 0.2 < x < 0.5 a 
ferromagnetic (FM) metallic behavior is usually found. 
As illustrated in Fig. [T] different magnetic orderings lead 
to diverse orbital configurations. The antiferromagnetic 
configurations are also insulating while fcrromagnetism is 
usually accompanied by metallicity by virtue of the DE 
interaction. 
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Manganite multilayers and thin films are epitaxially 
grown on insulating substrates. The lattice mismatch 
between the different layers gives rise to a uniaxial strain 
which may affect the bulk properties of a thin mangan- 
ite^ZHSZl and/or the electronic reconstruction at its inter- 
face .^ St rain may also be responsible for phase separa- 
tion! 33 * 3 ^ The inplane strain may range from ~ —2.3% to 
~ 3.2% depending on the substrate and the growth direc- 
tion.^ This strain modifies the relation between the lat- 
tice parameter in the direction perpendicular to growth 
(c) and the one in the parallel plane (o). 

The possible modifications produced by strain in a 
manganite are twofold!^ (i) the reduction (increase) of 
the lattice parameter in a particular direction would lead 
to an increase (reduction) of the hopping amplitude and, 
(ii) a distortion of the pseudocubic symmetry leads to a 
splitting of the e g levels which ma y produ ce orbital order- 
ing via the Jahn- Teller coupluig P 3 l 31 l 36 t BSl These mech- 
anisms alter the competition between the localizing and 
delocalizing interactions in manganites in opposite direc- 
tions. For instance, in a (001) thin film, a compressive 
strain (reduction of the lattice parameter in the xy plane) 
would increase the hopping within the xy plane and de- 
crease it in the z-direction favoring the x 2 — y 2 e g or- 

2 orbitals) to be occupied. 



bitals (rather than the 3z 
On the other hand, a compressive tetragonal distortion 
would produce a lowering o f the 3z 2 — r 2 orbitals with 
respect to the x 2 — y 2 ones! 23 1 31 ' 36 1 37 1 Experimentally,^ 
it is observed that when the lattice parameter in the z- 
direction is larger than in the xy plane (namely, c/a > 1) 
a C-type AF ordering with occupied 3z 2 — r 2 orbitals is 
favored in contrast to the A-type AF configuration with 
occupied x 2 — y 2 orbitals which occurs when c/a < 1. 
These observations are consistent with a predominance 
of mechanism (ii) over (i). 

Here we analyze the effect of strain on homogeneously 
strained epitaxially grown manganites in the (001) direc- 
tion. In Sec. |n]we study the modifications produced by 
strain in the phase diagram of half-doped (x = 0.5) bulk 
manganites. In Sec. |III| we turn the focus to the interface 
layers between a ferromagnetic metallic manganite with 
x = 0.3 and the insulating substrate and study how the 
electronic reconstruction is affected by strain. We fin- 



ish in Sec. IV with a discussion on the light of reported 
experimental results and the conclusions. 



II. EFFECT OF STRAIN IN BULK 
MANGANITES AT HALF-DOPING 



A. Model 

In the absence of strain, the model Hamiltonian for 
manganites includes the kinetic energy, the electron- 
lattice Jahn- Teller coupling, and the nearest neighbor 
antiferromagnetic superexchange^ 



+ E + + Qli + A (QuPi + Qht? + Q m t?)] 

i 

+ y.j^Sj (i) 



where C\ creates an electron on the Mn i-site in the e g 
orbital 7 (7 = 1,2 with 1 = |x 2 -y 2 ) and 2 = |3z 2 -r 2 )). 
The hopping amplitude fij depends on the Mn core spins 
orientation given by the spherical angles 9 and tf> via the 
double-exchange mechanism 

fi,j = cas(0</2) cos(0,-/2)-|-exp[z(V>i-%)] sin(0i/2) sin(6» J /2) , 

(2) 



and on the orbitals involved t 



*(y) 



+ x(y) 
,2 



22 — to' y > an d #2.2 = t„ where the su- 
perindices x,y, and z refer to the direction in the lattice. 
In the absence of strain t x Q ,v = t and f% — 4/3t. All the 
parameters are given in units of t which is estimated to 
be ~ 0.1 — 0.5 eV depending on the chemical compo- 
sition.^! Q H) with / = 1,2,3, are the phonon modes of 
the oxygen octahedra. pi — C\ is the site oc- 
cupation. The orbital occupation is represented by the 
Pauli matrices r I,z : rf — 2Re(Cj 3z2 _ r2 C ix 2_ y 2) and 

t* = Cl x2 _ y2 Ci^2_ y 2 - cl 3z2 _ r2 C.^ 3z 2_ r 2. f3 is assumed 
to be very large so the Qi breathing mode is frozen. A 
is the Jahn- Teller coupling and Jaf is the antiferromag- 
netic superexchange coupling between the localized t2 9 
spins which is estimated to be ~ 1 — 10 meV [T]. 

The strain is introduced uniformly in a cubic system 
(c/a = 1). We consider uniaxial strain arising from the 
lattice mismatch between the manganite and a cubic sub- 
strate and assume a (001) growth direction. The strain 
can be tensile (extension in the xy plane and compres- 
sion in the z direction: c/a < 1) or compressive (com- 
pression in the xy plane and extension in the z direction: 
c/a > 1). The strain in the xy-plane e xy is defined as 
e xy = { a s ~ a ) /o-s with a the average lattice parameter 
of the manganite and a s the in-plane lattice parameter 
of the substrate. Therefore, e xy > (< 0) corresponds 
to tensile (compressive) strain. The relation between the 
strain in the z-direction e z and e xy is given by the Poisson 
rati o v as e z = —Ave xy with 0.3 < v < 0.4 for mangan- 
ites | i2 | 30 | 37 | We choose e z = 

range between -0.02 (-2%) and 0.02 (2%) 

The effect of the strain on the system is twofold. On 
one hand, it affects the overlapping matrices and, there- 
fore the hopping amplitudes as-^1 



\e xy [41] and allow e xy to 



t z 



t(l - 2e xy ) 
4/3i(l-2e z ). 



(3) 



The strain can also induce a splitting S of the e g or- 
bitalsP323E3 tensile strain lowers the energy of the 
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FIG. 2: Bulk phase diagram A versus Jaf for half-doped manganites. FM, A, C and CE label the different magnetic and 
orbital orders considered (see Fig. [TJ. (a) Without strain, (b) With compressive strain e xy = —2%. (c) With tensile strain 
e xy — 2%. The splitting between the e 3 levels is 8 = 50e xy t [namely, \8\ — t in (b) and (c)]. 




FIG. 3: Bulk phase diagrams Jaf versus e xy for A = It with 8 = (a) and 8 — 100e xy t (b). 



\x 2 — y 2 ) orbital with respect to \Zz 2 — r 2 ) (t z > 0) while 
compressive strain does the opposite (t z < 0). This is 
introduced as an extra term in the Hamiltonian 



7 i 



with e 3z 2_ r 2 = 8/2 and e x -2_ y i = —5/2. We have an- 
alyzed a range of values \S\ < I00\e xy \t (namely \S\ < 
2t)P 

We find the ground state configuration of a half-doped 
manganite by solving the Hamiltonian in Eq. [l] plus the 
term in Eq. [4] self-consistently (at zero temperature) in 
a 4 x 4 x 4 system with periodic boundary conditions in 
the three directions. The phase diagrams as a function 
of A, Jaf an d e xy result from comparing the energies of 
the different configurations illustrated in Fig. [l] 



B. Results 

The phase diagram as a function of A and Jaf for a 
half-doped manganite is well knowrP^l and is shown in 
Fig. [2ja) as a reference. For small values of Jaf and 
A < 1.6£, the ground state is ferromagnetic and metallic. 
For A > 1.6t (not shown), the FM phase is insulating. As 
Jaf increases, and for A < 0.7t, two different antiferro- 
magnetic phases arise: (i) For a narrow range of values 
of Jaf, an A- type AF phase consisting of FM xy-planes 
coupled antifcrromagnetically in the z-direction, and (ii), 
for the largest realistic values of Jaf, the CE AF order. 
The CE order consists of FM zig-zag chains coupled AF 
between them. As A increases, this CE phase becomes 
dominant. 

In Fig. [3] the effect of the strain on the ground state 
configuration is shown for A = It. In Fig. [3|a), the 
splitting 8 between the e 9 levels caused by the strain 
is neglected so the only effect of the strain is to modify 
the hoppings as described in Eq. [3] . Compressive strain 
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(e xy < 0) produces a decrease of the hopping in the z- 
direction and an enhancement in the xy-plane causing 
the FM configuration to lose energy with respect to the 
A-AF configuration. On the contrary, tensile strain in- 
creases the hopping in the z-direction and lowers it in the 
xy-plane. In this case, both AF phases (A and CE), with 
zero hopping in the z-direction due to the antiferromag- 
netic order, lose kinetic energy with respect to the FM 
phase. 

When the splitting S is included (Eq. [4| , the changes 
in the phase diagram are more dramatic, see Fig. [3jb). 
As illustrated in Fig. [T] the different AF configurations 
are related to specific orbital ordcrings. In terms of the 
pseudospin r z , the C-AF ordering has t z < (preferred 
occupation of the 2>z 2 — r 2 orbital) while both the A and 
CE orderings have r z > (preferred occupation of the 
x 2 — y 2 orbital). The term in Eq. [4] tends to enforce a 
particular value of t z . Therefore, the A phase becomes 
the ground state for tensile strain (5 > 0) due to the 
lowering of the x 2 — y 2 orbital with respect to the 3z 2 — 
r and a C phase dominates for compressive strain and 
Jaf ^ O.K. This C phase appears for \5\ > 50\e xy \t 
if A = and for |<5| > 60|e.,; y |i if A = It. Comparing 
Figs. [3]ja) and (b), it is apparent that the effect of the 
strain on the splitting of the e g levels clearly overcomes 
the effect caused by the modifications in the hopping for 
a sufficiently large 5 > 50\e xy \t. 

In Fig. [2] (b) and (c) the modifications on the A versus 
Jaf phase diagram caused by compressive and tensile 2% 
strain are illustrated with \5\ = t. Consistently with the 
results in Fig.[3jb), the C phase arises when compressive 
strain is applied, while the A phase becomes more promi- 
nent with tensile strain. It is also noticeable in Fig. [2] (b) 
that the C ordering is not favored by the JT coupling A. 



III. EFFECT OF STRAIN AT A 
M ANG ANITE /INSULATOR INTERFACE 

A. Model 

We turn now to analyzing the effect of the strain on the 
electronic reconstruction at a (001) manganite/insulator 
interface. We focus on the case of a manganite with x — 
0.3, which corresponds to a ferromagnetic and metallic 
bulk ground state. In the (001) direction, manganitcs 
alternate MnC>2 and AO planes as shown in Fig. [4] The 
AO planes give a positive background charge 1 — x. We 
consider a thin manganite slab with a 4 x 4 cross section 
and l z — 12 Mn planes in the z-direction. The insulator 
is included as a hard-wall fixing the boundary condition 
at I = 1 and I = l z to zero charge density.^ The interface 
layers are I — 2 and I = l z — 1. 

The model Ha milto nian is the same as described for the 
bulk case in Sec. |II A| plus a Hartree term that takes into 
account the long range Coulomb interaction between all 
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FIG. 4: (Color online) 2-dimensional (xz plane) projection 
of the considered heterostructure. The circles represent the 
Mn sites while the squares are the A-sites (La, Sr, Ca, etc) 
shifted by (1, 1, 1)| with respect to the Mn. a is the lattice 
parameter. Full squares represent the Aq^A'q^O 2- plane, 
with a charge density of +0.7 per A atom. Empty squares 
correspond to the Sr 2+ 2 ~ planes where the charge density 
is zero. Periodic boundary conditions are applied in all three 
directions. 



the charges in the system? 3 * 16 * 19 * HHartrce takes the form 



it _ e 2 ^ / l (ni)( nj ) 1 ZjZj Zjjn,} \ 

(5) 

with Ri the position of the Mn ions, (m) — ^7(^,7^,7) 
the occupation number on the Mn i-site, eZi the charge of 
the A-cation located at Rf, and e the dielectric constant 
of the material. The relative strength of the Coulomb 
interaction is given by the parameter a = e 2 /aet ^ 1 — 
2 [TT], where e is the dielectric constant of the manganite. 

The charge density is ~ 0.7 in the central planes where 
the bulk values for the x = 0.3 manganite are recovered. 
However, close to the insulator, the charge density de- 
creases towards to fulfill the boundary condition. A re- 
distribution of charge occurs to screen the positive charge 
background and is controlled by the Coulomb parameter 
a. Experimentally, the charge transfer between different 
layers occurs within 2 — 3 unit cells, ^ where the charge 
density may be close to half-doping and the CE-type AF 
ordering. 

We consider different possible configurations at the 
manganite-insulator interface and compare their energies 
to define the phase diagrams as a function of A, Jaf, G xy , 
and aP31 (i) FM corresponds to all ferromagnetic planes; 
(ii) ICE stands for a configuration with all ferromagnetic 
planes except for a single CE plane at the manganite- 
insulator interfaces; and (iii) 2CE includes two CE planes 
at the manganite-insulator interfaces. 
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FIG. 6: (Color online) Phase diagram a versus Jaf for the 
manganite/insulator interface with A = 1. FM stands for all 
FM Mn planes, ICE stands for a configuration with a single 
CE layer at the manganite surface, and 2CE stands for two 
CE layers at the manganite surface. The lines represent the 
boundaries between the different ground state configurations: 
The black lines (squares) correspond to the results with e xy = 
0, the red lines (stars) are for compressive strain e xy — —2%, 
and the blue lines (circles) for tensile strain e xy = 2%. 5 = 
50e xy t is assumed. 



FIG. 5: (Color online) Average charge per plane for the three 
possible configurations considered at a manganite/insulator 
interface for A = 1, a = 1 and e xy = 0. The straight lines 
in black represent the positive background charge, (a) FM: 
all ferromagnetic planes, (b) ICE: all ferromagnetic planes 
except for a single CE plane at each manganite/insulator in- 
terface (layers I = 2, 11), and (c) 2CE: two CE planes at each 
manganite insulator interface (layers 1 — 2,3, 10, 11). 



B. Results 

In Fig. |5] the average charge per plane in an insula- 
tor /manganite/insulator trilayer is shown for the three 
possible ground state configurations at the interface for 
A = It, a = 1 and e xy = 0. The redistribution of charge 
depends only very weakly on the value of the strain (not 
shown). The CE order opens a gap in the density of 
states at x = 0.5 so, in the case of a = it would tend 
to pin the charge density to this value However, for 
a finite a, the Coulomb term Eq. [5] controls the charge 
redistribution in such a way that the charge density is 
not pinned at 0.5 in the CE layers and it is controlled in- 
stead by the screening of the positive background charge 
by means of the Hartree term in Eq. [5j 

Fig. [6] shows the a versus Jaf phase diagram for three 
different values of the strain: 0, 2% and —2% with A = It 
and S — 50e xy t. This phase diagram without strain 
was studied before in Ref. |13) with a model which ne- 
glected Jahn- Teller interactions but considered instead 
an interorbital Hubbard term to stabilize the antiferro- 
magnetic phases. Those results compare very well to 




FIG. 7: (Color online) Phase diagram A versus Jaf for the 
manganite/insulator interface with a = 1. The labels are the 
same as in Fig. [6] 8 = 50e xy t is assumed. 

our e xy — results in Fig. [6] (squares). A compressive 
(e xy < 0) strain makes the FM configuration relatively 
more stable (namely, a larger Jaf is required to produce 
an antiferromagnetic CE order in one or two layers) while 
a tensile (e xy > 0) strain lowers the energy of the antifer- 
romagnetic CE interface configurations with respect to 
the FM. These results are consistent with the behavior 
observed for bulk manganites in Figs.[2]and[3l(b) and are 
due to the preferred occupation of the x 2 — jr orbitals in 
the CE order. 

Similar phenomenology is observed in Fig. [7] where 
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FIG. 8: Phase diagrams Jaf versus e xy for the manganite/insulator interface with A = 1 and a = 1. The strain induced e g 
level splitting is S — in (a) and 8 — 100e xy t in (b). 



the A versus Jaf phase diagram is plotted. In general 
terms, tensile strain tends to favor the antiferromagnetic 
CE order close to the insulator while compressive strain 
does not seem to affect the interface ground state very 
strongly. Finally, we show the Jaf versus e xy phase dia- 
gram in Fig. [8] for A = It and a = 1. In Fig. [8] (a) 8 — 
so the only effect of the strain is to modify the hopping 
amplitude as given in Eq. [3] The effect of the strain 
is very mild in this case with only a slight gain of the 
configurations with CE planes for compressive (e xy < 0) 
strain. The tendency is the opposite and the dependence 
on strain is stronger when the splitting 8 of the e g levels 
is included, as illustrated in Fig. [8jb) for 8 = I00e xy t. 

At interfa ces orbital reconstruction may 
arise JHH ?5|44|45 | j n particular, due to the breaking 
of the translational symmetry in the z-direction, the 
3z 2 — r 2 orbital cannot gain much kinetic energy 
producing a splitting of the e g levels which favors the 
occupation of the x 2 — y 2 orbital. In the all FM layers 
case, in the absence of strain, we observe this kind of 
ferro-orbital configuration (which corresponds to r 2 > 0) 
at the interface layer. A positive strain (e xy > 0) 
enhances this ferro-orbital interface ordering which also 
occurs in the bulk, as discussed in Sec. |TT] On the other 
hand, a negative strain would produce an e g splitting 
opposite to the preferred one at an interface, leading 
for a sufficiently large value of the strain to a reduction 
and, eventually, to a sign change of t 2 . Therefore, r 2 
or, equivalently, the Q3 phonon mode, increases as the 
strain goes from compressive to tensile. t x is nonzero 
only for the CE phase and has opposite signs on the 
sites with 3x 2 — r 2 and 3y 2 — r 2 -like orbitals (see Fig. [T|. 
\t x \ is a measure of the mixing of the two e g levels and 
therefore decreases as t z increases. 

CE layers adjacent to a FM layer, as we have in both 
the ICE and the 2CE configurations, show a charge and 
orbital distribution which is different from the isolated 
case. In particular, half of the spins are parallel to the 
spins in the nearby FM layer which produces a larger 



occupation on those sites (both on the CE and the FM 
layers) with respect to the sites with antiparallel spinsP^ 
Moreover, the redistribution of charge at the interface 
controlled by the Hartree term leads, in general, to a 
charge density away from the 0.5 which stabilizes CE 
(see Fig. [5]). As a consequence, in our calculations the CE 
phase, which is insulating in bulk at x — 0.5, may become 
metallic. We have analyzed the possibility that insulating 
behavior may arise at interfacial layers due to strain. We 
found that the interface CE layer becomes insulating (a 
gap opens at the Fermi energy) for relatively large values 
of tensile strain e xy ~ 2% and 8 = 2t if a < 0.7. For 
larger values of a and/or smaller or negative values of 
the strain, the interface CE layer is always metallic. In 
reality, we expect that phase separation^ or a different 
(maybe incommensurate with the lattice) ordei^ may 
arise at the interface layers leading to insulating behavior. 
In this case, our results imply that an already existing 
insulating gap may be enhanced by tensile strain. On 
the other hand, compressive strain favors the occupation 
of the 3z 2 — r 2 orbitals which gain energy by hopping 
in the z-direction to the adjacent FM layer, leading to 
metallic behavior. 



IV. DISCUSSION AND CONCLUSIONS 

It is well known that manganitc thin films have prop- 
erties (magnetic critical temperature T c , transport) dif- 
ferent from those exhibited by bulk manganites. In par- 
ticular, it has been observed that the conductivity and 
the magnetic T c are reduced with respect to bulkPSHU 
These differences may come about due to the effect of the 
strain, produced by the lattice mismatch with the sub- 
strate, throughout the whole film, and due to the modifi- 
cations at the substrate/manganite interface (electronic 
reconstruction, phase separation) , which can be most im- 
portant for the very thin films used in multilayers. The 
electronic reconstruction may also be strongly affected by 
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strain as it is mainly related to a redistribution of charge 
(although other effects, like disorder, might also play a 
role) . We have focused here on the effect of the strain on 
the electronic properties of manganites at the interface 
with an insulator and in bulk. 

The "active" orbitals in manganites are the Mn e s : 



3z 2 



r 2 and x 2 



y . A tetragonal distortion of the 



(pseudo)cubic unit cell produces a preferred occupation 
of one of the anisotropic e g levels or, equivalently, to 
an energy splitting which is associated to a particular 
Ja hn- Teller p honon mode and leads to orbital order- 
m g) 2 3 | 3i | 36 | 37 | j n turn, a particular orbital ordering is as- 
sociated to a particular magnetic ordering (see Fig. [I]). 
Therefore, an extension of the lattice parameter in the 
xy-plane (tensile strain) favors the occupation of x 2 — y 2 
which is dominant in the A and CE-type AF orders while 
compressive strain favors 3z 2 — r 2 and, hence, the C-type 
order. This is exactly what we find for x = 0.5, see Figs. [2] 
and [3j if the splitting 5 is sufficiently large (S > 50e xy t) . 
This value of S is relatively large (~ t for 2% strain) but 
comparable to the Jahn- Teller splitting.^ This tuning of 
the orbital arran gement with strain has been found in 
experiments) 12 * 36 ^ model calculations for LaMn03p2l and 
ab-initio calculations for Lao.66Sro.33MnC>3 jH] and for 
LaMn0 3 /SrMn0 3 superlatticesP 

The general trends are the same in the insula- 
tor/manganite (x — 0.7)/insulator sandwich considered 
in Sec. Ill The bulk behavior (FM and metallic) is recov- 
ered within a few unit cells while at the interface layer, 
where a charge density close to 0.5 is expected due to 
the redistribution of charge, a CE order may arise. This 



CE order is favored by tensile strain and disfavored by 
compressive strain, as in bulk. 

Due to double exchange, FM and metallicity usually 
come hand in hand in manganites. DE is suppressed 
when the degeneracy of the e ff levels is broken. Therefore, 
it is expected that uniform st rain produces a reduction of 
the T C .E3 As discussed in Sec. Ill at an interface layer the 



splitting of the e g levels produced by compressive strain 
(S > 0) is enhanced. This may also lead to a reduced 
conductance at interfaces with respect to bulk. 

We are assuming a uniform strain and a uniform Pois- 
son ratio v. However, even in the case of achieving a 
uniform in-plane strain, the strain in the z-direction may 
change with the distance to an interface or surface. In 
this case, 8 may have different signs at different atomic 
planes and lead to more complex orbital arrangements 
than the ones reported here.^ 

In conclusion, we have performed model calculations of 
the effect of uniaxial uniform strain on a bulk half-doped 
manganite and at the interface between a ferromagnetic 
metallic manganite and an insulator. The main conse- 
quence of the strain is the occurrence of an orbital se- 
lection which is intimately related to the spin degree of 
freedom and the transport properties. In this way, strain 
provides with a way to tune the ground state configu- 
ration on manganites and, therefore, control the perfor- 
mance of manganite based electronic devices. 

This work is supported by FIS2009-08744 (MICINN, 
Spain). A.B. also acknowledges the JAE program 
(CSIC, Spain) and M.J.C. the Ramon y Cajal program 
(MICINN, Spain). 



1 E. Dagotto, Nanoscale Phase Separation and Colossal 
Magnetoresistance (Springer- Verlag, Berlin, 2003). 

2 C. Israel, M. J. Calderon, and N. D. Mathur, Materials 
Today 10, 24 (2007). 

3 M.-H. Jo, N. D. Mathur, N. Todd, and M. G. Blamire, 
Phys. Rev. B 61, R14905 (2000). 

4 M. Bibes and A. Barthelemy, IEEE Trans. Electron. De- 
vices 54, 1003 (2007). 

5 E. Dagotto, Science 318, 1076 (2007). 

6 M. Huijben, A. Brinkman, G. Koster, G. Rijnders, 
H. Hilgenkamp, and D. H. A. Blank, Adv. Mater. 21, 1665 
(2009). 

7 A. Ohtomo, D. Muller, J. Grazul, and H. Hwang, Nature 
419, 378 (2002). 

8 A. Ohtomo and H. Y. Hwang, Nature 427, 423 (2004). 

9 S. Okamoto and A. Millis, Nature 428, 630 (2004). 

10 T. Koida, M. Lippmaa, T. Fukumura, K. Itaka, Y. Mat- 
sumoto, M. Kawasaki, and H. Koinuma, Phys. Rev. B 66, 
144418 (2002). 

11 C. Lin, S. Okamoto, and A. J. Millis, Phys. Rev. B 73, 
041104 (2006). 

12 H. Yamada, M. Kawasaki, T. Lottermoser, T. Arima, and 
Y. Tokura, Appl. Phys. Lett. 89, 052506 (2006). 

13 L. Brey, Phys. Rev. B 75, 104423 (2007). 

14 §erban Smadici, P. Abbamonte, A. Bhattacharya, X. Zhai, 
B. Jiang, A. Rusydi, J. N. Eckstein, S. D. Bader, and J.-M. 



Zuo, Physical Review Letters 99, 196404 (2007). 

15 D. Niebieskikwiat, M. B. Salomon, L. E. Hueso, N. D. 
Mathur, and J. A. Borchers, Phys. Rev. Lett. 99, 247207 
(2007). 

16 J. Salafranca, M. J. Calderon, and L. Brey, Phys. Rev. B 

77, 014441 (2008). 

17 S. J. May, A. Shah, S. G. E. te Velthuis, M. R. Fitzsim- 
mons, J. M. Zuo, X. Zhai, J. N. Eckstein, S. D. Bader, and 
A. Bhattacharya, Phys. Rev. B 77, 174409 (2008). 

18 C. Lin and A. J. Millis, Phys. Rev. B 78, 184405 (2008). 

19 M. J. Calderon, J. Salafranca, and L. Brey, Phys. Rev. B 

78, 024415 (2008). 

20 B. R. K. Nanda and S. Satpathy, Phys. Rev. Lett. 101, 
127201 (2008). 

21 R. Yu, S. Yunoki, S. Dong, and E. Dagotto, Phys. Rev. B 
80, 125115 (2009). 

22 Z. Sefrioui, C. Visani, M. Calderon, K. March, 
C. Carretero, M. Walls, A. Rivera-Calzada, C. Leon, 
R. Lopez Anton, T. R. Charlton, et al., Advanced Ma- 
terials (2010). 

23 L. Abad, V. Laukhin, S. Valencia, A. Gaup, W. Gudat, 
L. Balcells, and B. Martinez, Advanced Functional Mate- 
rials 17 (2007). 

24 J. Chakhalian, J. W. Freeland, H.-U. Habermeier, G. Cris- 
tiani, G. Khaliullin, M. van Veenendaal, and B. Keimer, 
Science 318 (2007). 



8 



25 P. Yu, J.-S. Lee, S. Okamoto, M. D. Rossell, M. Huijben, 
C.-H. Yang, Q. He, J. X. Zhang, S. Y. Yang, M. J. Lee, 
et al., Phys. Rev. Lett. 105, 027201 (2010). 

26 R. Kajimoto, H. Yoshizawa, Y. Tomioka, and Y. Tokura, 
Phys. Rev. B 66, 180402 (2002). 

27 Z. F. et al., Phys. Rev. Lett. 84, 3169 (2000). 

28 Y. Ogimoto, M. Izumi, T. Manako, T. Kimura, 
Y. Tomioka, M. Kawasaki, and Y. Tokura, Appl. Phys. 
Lett. 78, 3505 (2001). 

29 I. Infante, S. Estrade, F. Sanchez, J. Arbiol, F. Peiro, 
V. Laukhin, J. P. Espinos, M. Wojcik, E. Jedryka, and 
J. Fontcuberta, Phys. Rev. B 76, 224415 (2007). 

30 C. Adamo, X. Ke, H. Q. Wang, H. L. Xin, T. Heeg, M. E. 
Hawley, W. Zander, J. Schubert, P. Schiffer, D. A. Muller, 
et al., Appl. Phys. Lett. 95, 112504 (2009). 

31 A. Sadoc, B. Mercey, C. Simon, D. Grebille, W. Prellier, 
and M.-B. Lepetit, Phys. Rev. Lett. 104, 046804 (2010). 

32 B. R. K. Nanda and S. Satpathy (2009), arXiv:0911.0683. 

33 K. H. Ahn, T. Lookman, and A. R. Bishop, Nature 428, 
401 (2004). 

34 P. R. Sagdeo, R. J. Choudhary, and D. M. Phase, J. Appl. 
Phys. 104 (2008). 

35 A. J. Millis, T. Darling, and A. Migliori, J. Appl. Phys. 
83, 1588 (1998). 

36 Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 



37 B. R. K. Nanda and S. Satpathy, Phys. Rev. B 78, 054427 
(2008). 

38 S. Dong, S. Yunoki, X. Zhang, C. §en, J.-M. Liu, and 
E. Dagotto, Phys. Rev. B 82, 035118 (2010). 

39 L. Brey, Phys. Rev. Lett. 92, 127202 (2004). 

40 A. Ramirez, J. Phys.: Condens. Matter 9, 8171 (1997). 

41 K. H. Ahn and A. J. Millis, Phys. Rev. B 64, 115103 
(2001). 

42 W. Harrison, Electronic Structure and the Properties of 
Solids: The Physics of the Chemical Bond (Dover Publi- 
cations Inc, New York, 1989). 

43 L. Brey, Phys. Rev. B 71, 174426 (2005). 

44 M. J. Calderon, L. Brey, and F. Guinea, Phys. Rev. B 60, 
6698 (1999). 

45 S. Okamoto, Phys. Rev. B 82, 024427 (2010). 

46 G. Milward, M. J. Calderon, and P. B. Littlewood, Nature 
(London) 433, 607 (2005). 

47 I. C. Infante, F. Sanchez, J. Fontcuberta, M. Wojcik, 
E. Jedryka, S. Estrade, F. Peiro, J. Arbiol, V. Laukhin, 
and J. P. Espinos, Phys. Rev. B 76, 224415 (2007). 

48 A. Millis, Phys. Rev. B 53, 8434 (1996). 

49 C. Ma, Z. Yang, and S. Picozzi, J. Phys.: Condens. Matter 
18 (2006). 



